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We consider site percolation on a correlated bi-colored simple cubic lattice. The corre- 
lated medium is constructed from a strongly alternating bi-colored simple cubic lattice 
due to anti-site disordering. The percolation threshold is estimated. The cluster size dis- 
tribution is obtained. A possible application to the double 1:1 perovskites is discussed. 
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1. Introduction 

In percolation theory one aims to model random (disordered) media to evaluate their 
typical properties. 1 The simplest example is called Bernoulli percolation, where all 
sites have the same probability to be occupied independent of each other. The 
assumption of stochastic independence is too strong for consider many natural sys- 
tems. For such systems a medium is not Bernoulli rather then correlated. Recently, 
a particular lattice topology and a particular constraint were investigated. 2 ' 3 We 
will deal with another particular correlated medium. 

Let us consider a strongly alternating bi-colored simple cubic lattice (Fig. 1). 
Each 'white' site is surrounded by 'black' ones and vice versa, i.e. each of the sites 
is isolated and forms a cluster of unit size. The medium is fully ordered. One 
can divide the lattice into two sublattices: 'black' and 'white'. The probability 
that a 'white' site belongs to the 'white' sublattice is unity. Let us assume that 
alternating order is partially destroyed. We will denote p the probability to find a 
'white' site in the 'black' sublattice. If p > 0, then one can find a 'white' cluster 
of size more than the unity. If p = 0.5 (fully disordered system) then there is an 
infinity cluster indeed, because the site percolation threshold for the simple cubic 
lattice and Bernoulli percolation is 0.31 1608. 4 We will look for a critical disturbance 
p c where an infinity cluster is formed, for the first time. 



2. Cluster size distribution 
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Fig. 1. Bi-colorcd alternating simple cubic lattice (rock salt lattice). 

The probabilities per cluster site n s that a randomly chosen lattice site belongs to 
a cluster of s sites are 

n 1 = \{{l-pf+p 7 ), (1) 
n 2 = 3/(l-/ 6 , (2) 
n 3 = |p 5 (l - pf (p 7 + (1 - p) 7 + 4p* + 4(1 - pf) , (3) 

n i = 3p 11 (l-p) 11 +51p 10 {l-p) 10 

+ 6/(1 - pf (p 12 + (1 - p) 12 ) + 4/(1 - pf (p 11 + (1 - p) 11 ) , (4) 

n 5 - \p w (l - p) w ( P 7 + (1 - p) 7 + 8/ + 8(1 - pf 

+32/ + 32(1 - pf + 8/ + 8(1 - pf) 
+ 12/(1 - pf (/ + (1 - Z 8 + 3/ + 3(1 - pf + 8 (/ + (1 - pf)) 

+ y/(l-p) 3 (f 17 + (!-f) 17 )- (5) 

The numerical results were obtained using Hoshen-Kopelman algorithm. 5 We 
investigated a number of sample lattices with linear size L up to 128 sites. The 
theoretical results and the computer simulations are in the reasonable agreement 
(see Fig. 2-6). 

Cluster size distribution is surprising. In contrast with Bernoulli percolation 
one can see the oscillations (see Fig. 7-8). If p = 0.5 (disordered medium) then the 
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distribution demonstrates classical behavior (Fig. 9). In Section 4 we will discuss 
a possible connection such oscillations with the experiments and the real world 
systems. 
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Fig. 2. The probability that an arbitrary site belongs to cluster of size s = 1. Solid line — from 
Eq. 1, crosses — numerical data. Here and below only left part of the plot is presented because of 
symmetry. 



3. Numerical estimation of percolation probabilities 

Estimates for p c have been obtained by means of percolation frequencies. Simu- 
lations give the percolation frequencies, which serve as an approximation of the 
percolation probability. Critical percolation have been estimated by not-linear fit 
functions defined by 

p=i(l + erf(( P -p c )/a)) (6) 

and 

p = (l + exp(-(p-p c )a)) _1 . (7) 

The first function converges to a step function as a — * and the second one as 
a — ► oo. Moreover we used a polynomial of 3-th degree. The fits gave almost the 
same results for the places (the values of p) where the spanning probability curve 
reaches a value of 1/2. The similar approach to estimate the critical probabilities 
was performed in works. 2,6 Numerical estimation suggests that p c « 0.147 for L = 
128 (Fig. 10). 

The percolation threshold was calculated for three values of the linear lattice 
size L — 32, 64, 128. The percolation threshold p c (oo) for infinite lattices can be 
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Fig. 3. The probability that an arbitrary site belongs to cluster of size s = 2. Solid line — 
from Eq. 2, crosses — numerical data. Here and below one can see a visible deviation from the 
theoretical curve because of well known finite size effect for the large clusters. 1 
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Fig. 4. The probability that an arbitrary site belongs to cluster of size s = 3. Solid line — from 
Eq. 3, crosses — numerical data. 
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Fig. 6. The probability that an arbitrary site belongs to cluster of size s 
Eq. 5, crosses — numerical data. 
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Fig. 7. Cluster size distribution for p = 0.05. Here and below cluster of unit size is omitted because 
of very large magnitude. 
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Fig. 8. Cluster size distribution for p = 0.1465. 
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Fig. 10. Percolation frequency as a function of p (L = 128) (high resolution). 
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found by fitting these results for different lattice sizes to the scaling relation 1 

|p c (L)-p c (oo)| «r^, (8) 

where the critical exponent v has the value 0.875 in three dimensions. 7 This method 
leads to an estimate p c (oo) « 0.1454 (Fig. 11). 
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Fig. 11. The percolation threshold as a function of L x l v . 

We observed some infinity clusters at percolation threshold, nevertheless we 
cannot assert there is not a finite size effect. 8 

4. Possible applications to the double 1:1 perovskites 

In last decades, much attention has been paid to ABO3 perovskite-type oxides due to 
their unique dielectric, electrooptic, and other properties. Complex perovskite-like 
oxides of the AB^B^^Os type, in particularly double 1:1 perovskites AB]y 2 B^ , 2 03* 
e.g. (Ba, Sr)Fe 1 / 2 Mo 1 / 2 3 and PbFei/ 2 Nbi/ 2 03, have much more surprising and 
useful properties. Examples of their applications are piezoelectric transducers and 
actuators, non-volatile ferroelectric memories, and microelectronic devices. The 
structure is built up by ordering perovskite blocks in a rock salt superlattice (Fig. 12) 
and the properties of the material are thought to critically depend on this ordering. 
The double perovskite Sr 2 FeMo06 and related materials are good candidates for 
magnetic devices, as they combine a high Curie temperature and a fully polarized 
(half metallic) conduction band. For more detailed overview and a large collection 
of references see e.g. 9 ' 10 . 



*Such the alloys are denoted as A2BB'C>6, too. 
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Fig. 12. An example of ordered 1:1 double perovskite: PbFe^Nb^Oa. Atoms of Fe and Nb 
form an alternating simple cubic lattice. 

Little thing is known and understood about the microscopic mechanisms re- 
sponsible for their very convenient properties. At the present, these materials are 
being extensively studied. One particular aspect that seems very promising is the 
effect of atomic ordering on the properties of these alloys. Anti-site disordering in 
Sr 2 FeMo06 double perovskites (containing B atoms at B' positions, and vice versa) 
has recently been shown to have a dramatic influence in their magnetic and magne- 
totransport properties. Experiments suggest that, in many samples, the saturation 
magnetization is less than the expected value. This effect is usually ascribed to the 
presence of antisite defects, where, due to the similarity of their atomic radii, Mo 
ions are randomly placed on the Fe sublattice and conversely. 

A model described above seems very promising to understand some features 
of the double perovskites. In particularly, the oscillations of the hypcrfine mag- 
netic fields distribution in the PbFei^Nb^Oa samples 11 may be related with the 
oscillations of the cluster size distribution (Fig. 7,8). 

Acknowledgements 

The authors are grateful to Prof. R. V. Vcdrinskii for pointing out the problem. 
They thank Prof. S. A. Prosandeev and Prof. I. P. Raevskii for discussions. 

References 

1. D. Stauffer, A. Aharony, Introduction to Percolation Theory (Taylor & Francis, London, 
1992). 

2. S. Reimann, J. Bendisch, Physica A 296, 391 (2001). 

3. J. Bendisch, S. Reimann, H. von Trotha, Physica A 307, 1 (2002). 



Yu. Yu. Tarasevich & E.N. Manzhosova 



4. C. D. Lorenz, R. M. Ziff, J. Phys. A: Math. Gen. 31, 8147 (1998). 

5. J. Hoshen, R. Kopelman, Phys. Rev. B 14, 3438 (1976). 

6. F. Babalievski, Phys. Rev. E 55, 1228 (1997). 

7. A. Bunde, S. Havlin in Fractals and Disordered Systems ed. A. Bunde, S. Havlin 
(Springer, 1996), p. 65. 

8. D. Stauffer, Physica A 242, 1 (1997). 

9. L. Balcells, J. Navarro, M. Bibes, et al., cond-mat/0007248. 

10. J. L. Alonso, L. A. Fernandez, F. Guinea, et al., cond-mat/02 10303. 

11. I. P. Raevskii, V. V. Kitaev, S. A. Brugerman, et al., in Order, Disorder and Properties 
of the Oxides — Proc. Sochi Int. Symp. (2002, v. 2), p. 54. 



